Surface-enhanced pair transfer in quadrupole states of neutron-rich Sn isotopes 
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We investigate the neutron pair transfer modes associated with the low-lying quadrupole states 
in neutron-rich Sn isotopes by means of the quasiparticle random phase approximation based on 
the Skyrme-Hartree-Fock-Bogoliubov mean field model. The transition strength of the quadrupole 
pair-addition mode feeding the 2^ state is enhanced in the Sn isotopes with A > 132. The transition 
density of the pair-addition mode has a large spatial extension in the exterior of nucleus, reaching 
far to r ~ 12 — 13 fm. The quadrupole pair-addition mode refiects sensitively a possible increase of 
the effective pairing interaction strength in the surface and exterior regions of neutron-rich nuclei. 
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I. INTRODUCTION 

The pair correlation is one of the fundamental nucleon 
many-body correlations, which governs many aspects of 
the nuclear structure [l|-|l|. Recently attentions are paid 
to the properties of the pair correlation in neutron-rich 
nuclei, which exhibit new features such as the neutron 
halo, the neutron skin and the presence of weakly bound 
neutrons. For example, a spatial two-neutron correlation 
or the di-neutron correlation, between two weakly-bound 
neutrons forming the halo in ^^Li and ^He has been 
discussed extensively [S^dOij. and its experimental signa- 
tures are reported recently [ll|, [l2|. There exists also 
several theoretical studies demonstrating similar spatial 
two-neutron correlation in many other nuclear systems. 
The spatial two-particle correlation was originally dis- 
cussed in the closed-shell + two-particle systems [l3|-[l6{ , 
and recent selfconsistent Hartree-Fock-Bogoliubov mean- 
field calculations suggest that the correlation is gener- 
ically enhanced around the nuclear surface of neutron- 
rich nuclei covering also medium- and heavy-mass regions 
|17h19| . Furthermore it is found that neutron matter at 
low densities may exhibit features of the strong coupling 
pairing |20l - t22| . characterized by the small size (~ 5 fm at 
minimum) of the Cooper pair comparable with the av- 
erage inter- neutron distance 23-25], refiecting the large 
scattering length a = — 18.5fm of the bare nuclear force. 
This may also suggest a possible enhancement of the pair 
correlation in the surface and exterior regions of neutron- 
rich nuclei. The spatial two-neutron correlation and the 
surface enhancement of the neutron pairing is also dis- 
cussed in the cases of a slab and a semi-infinite nuclear 
matter [26J, |27[. Also the pairing interaction induced by 
the coupling to the phonons is expected to increase the 
surface enhancement |28| . 

It is important to ask how one can probe the surface 
enhancement of the pair correlation in neutron-rich nu- 
clei. We consider here the transfer of a neutron pair 
since the matrix element of the pair transfer has direct 
relation to the degree of the pair correlation [1, 0, [23- 
ISlj . Recently the neutron pair transfer process in light 
neutron-rich nuclei ^^Li and ^'*He has been studied ex- 



perimentally |32H3q . In the theoretical side, a possible 
link between the pair transfer between the 0"*" ground 
states and the surface enhancement of the pairing in 
medium and heavy neutron-rich nuclei was suggested by 
Dobaczewski et al.\3&i. Recently the pair transfer feed- 
ing the excited 0+ states in neutron-rich nuclei has been 
investigated |37H39j from the viewpoint of the pair vi- 
brational mode, which are collective vibrational modes 
associated with the pair correlation [l[ [2, |30|. Khan et 
al. 39] has argued that the pair transfer can be a possible 
probe of different models of the pairing interactions with 
differences in density dependence and surface enhance- 
ment. 

In the present work, we shed light on the neutron pair 
transfer with the quadrupole multipolaritv 13 ll |40|. l4l| .i.e. . 
the transitions from the ground state of a nucleus with 
neutron number N feeding excited 2+ states, in particu- 
lar the first 2+ state, in the neighboring A^ ± 2 isotopes. 
It is noted that the first 2+ state in open shell nuclei 
is a collective state of the surface vibrational type with 
a strong infiuence by the pair correlation[l|. With this 
surface character of the mode, we expect that the pair 
transfer feeding the 2^ state may provide information on 
the pair correlation possibly enhanced in the surface and 
exterior regions of neutron-rich nuclei. 

The low-lying quadrupole vibration is revealed recently 
to show new features in some neutron-rich nuclei. For 
instance, the ratio between the neutron and proton am- 
plitude becomes significantly different from the liquid 
drop ratio iV/Z [42-113] ■ New aspects of the quadrupole 
vibration may emerge also in the pair transfer transi- 
tions feeding the 2+ states. In our previous studies 
on the dipole and octupole responses, we have shown 
that soft modes accompanying a large pair transition 
amplitude in the surface and exterior regions emerge 
above the threshold ener gy i n nuclei with small neu- 
tron separation energy [13, El, H^. We also found that 
these soft dipole and octupole modes are sensitive to the 
surface-enhancement of the pair correlation. The low- 
lying quadrupole vibration is different from these soft 
modes in many respects: it exists systematically from 
stable nuclei to nuclei far from the stability line, and the 



excitation energy is located below the separation energy 
in most of the cases. Nevertheless, as we shall show be- 
low, the low-lying quadrupole mode also has sensitivity 
to the surface pair correlation. The sensitivity is sys- 
tematically seen in the pair-addition transfer feeding the 
lowest-lying quadrupole state, and it becomes significant 
in neutron-rich nuclei beyond the shell closure, i.e. in Sn 
isotopes with A > 132. 

We shall perform systematic numerical study by adopt- 
ing the Sn isotopes as a representative case since its long 
isotopic chain enables us to see systematic trends. In 
particular we focus on isotopes beyond the TV = 82 shell 
gap, above which the neutron separation energy sud- 
denly decreases. Our study is based on the Skyrme- 
Hartree-Fock-Bogoliubov mean-field theory in describ- 
ing the pair correlated ground state, and the quasipar- 
ticle random phase approximation fQRPA)p8| (and ref- 
erences in [531) to describe the 2+ states and the associ- 
ated pair transfer modes. In the present work, we con- 
centrate on the nuclear structure aspects, and leave the 
cross section and the reaction mechanism of the neutron 
pair transfer [3l|, [5ll-l53| for future works. The formal- 
ism is described in Section 2. To elucidate effects of the 
surface-enhanced pairing, we compare different models 
of the effective pairing interaction, which have and do 
not have the property of the surface enhancement. After 
presenting the adopted effective pairing interactions in 
Section 3, we discuss numerical results in Section 4, and 
the conclusions are drawn in Section 5. 



II. SKYRME-HARTREE-FOCK-BOGOLIUBOV 
MEAN-FIELD PLUS QRPA APPROACH 

A. model 

We adopt the Hartree-Fock-Bogoliubov mean-field 
model and the quasiparticle random phase approxima- 
tion (QRPA) to describe the pair transfer modes with 
the quadrupole multipolarity. The basic framework 
of the present HFB+QRPA formalism is described in 
Refs.[48l. |54|. Here we recapitulate it briefly in order to 
explain the procedure to describe the pair transfer modes. 

In performing the HFB calculation, we assume the 
spherical symmetry for the ground state and the asso- 
ciated mean-fields. The Skyrme energy density func- 
tional (the energy density functional associated with the 
Skyrme interaction) is adopted to derive the Hartree- 
Fock potential and the position-dependent effective 
mass|55j. The parameter set SLy4 is adopted |56{. As 
the residual interaction responsible for the pairing cor- 
relation we use the density-dependent delta interaction 
(DDDI) |,i|2l,|3i,|I^i| 

vf'ir,r') = ^{l-P.)V,{r)6{r-r') {q = n,p) (1) 

where the the pairing interaction strength Vq{r) is a func- 
tion of the neutron and proton densities Pq{r) {q = n,p), 



and hence it depends on the position r. We solve the 
Hartree-Fock-Bogoliubov equation in the radial coordi- 
nate space to obtain the wave functions of the quasi- 
particle states. We assume the box boundary condition 
at r = i?max with i?max — 20 fm. We truncate the 
quasiparticle states by the maximum quasiparticle en- 
ergy E'max = 60 MeV. We also truncate the quasiparticle 
states by the maximum orbital angular quantum num- 
ber /,„ax = 12, which is suflficiently large to obtain the 
converged HFB solutions of the Sn isotopes. 

The QRPA calculation is then performed on the ba- 
sis of the self-consistent HFB mean-fields. As the resid- 
ual interaction in the particle-particle channel, we adopt 
the DDDI, Eq. ([T]), the same effective pairing interac- 
tion used for the HFB calculation. Concerning the resid- 
ual interaction in the particle-hole channel, we derive it 
also from the Skyrme functional used for the HFB cal- 
culation, but we apply the Landau-Migdal approxima- 
tion to the residual interaction[37, 39, 48, 50]. Namely, 
we approximate the particle-hole residual interaction by 
a position-dependent contact interaction Vph{r,r') — 
{{Fo/No){r) + {F^/No){r)T ■ r')} S{r - r') where Fo/No 
and Fq/Nq are the Landau-Migdal parameters derived 
from the Skyrme functional. The position dependence is 
taken into account in the Landau parameters by the lo- 
cal density approximation to the Fermi momentum. As 
a prescription commonly adopted in the QRPA studies 
using the Landau-Migdal approximation [33, |3^, Ea l50| . 
we introduce a renormalization factor / multiplying Vph, 
and it is fixed to impose the zero excitation energy for 
the lowest-energy isoscalar dipole solution corresponding 
to the spurious center-of-mass motion. 

We adopt the linear response formalism for the QRPA 
calculation 54| . Given the above residual interactions, 
we have coupled linear response equations for three 
kinds of fluctuating densities Spq{r,uj), Spq_+{r,uj) and 
5pq^-{r,io), which correspond respectively to fluctuating 
parts of the nucleon density 



Pqir,t) ^ {^it)\pqir)\^it)) 

= {^{t)\J2^l{ra)ijq{ra)\^{t)) 

and of two kinds of the pair densities 



(2) 



p,,±(r,i) = {^{t)\^l{r ;)Vt(r t)±V'g(r- t)V',(r ;)|*(t)). 

(3) 
The time-dependence of the state vector, and hence the 
density fluctuations 5pq{r,uj) = YLMif')^Pq{r,uj)/r^ and 
Spq^±{r,uj) = Ylm i'f')^ Pq,±ir, uj) /r'^ in the frequency do- 
main arise from a response of the system caused by ex- 
ternal flelds with multipolarity L. It is convenient to 
introduce operators 



Pq{r) 



2 

-^tljqira)il^q{ra). 



(4) 
(5) 



Here PUr) and Pq{r) represent respectively addition and 
removal at position r of a pair of like-nucleons whose 
spins are coupled to the singlet 5* = 0. The pair densi- 
ties are written compactly as Pq.±{r,t) = {'ii{t)\P^{r) ± 



Pq{r)\^it)). 

For a response with the multipolarity L and the natural 
parity tt — {—)^, the coupled linear response equation is 
written in the radial coordinate system as 



SpgL{r,uj) 
5p+,qL{r,uj) 
5p-,qLir,u}) 



dr'\ RZLir,r',u^) 



Eq' '^11 ir')Spq.Ur',u;)/r'' + z^-^Jr') 



vTirir') 



K^^i,{r')6p+^qL{r',uj)/r''- t u^^^y 
-Kp,Ur')6p.,qUr',u:)/r'^ + f}-'^Jr') 



(6) 



where Rq qh ^^ ^^^ unperturbed response function rep- 
resenting non-interacting two-quasiparticle excitations. 
KphiT) and Kpair('') are the residual interactions in the 
particle-hole and the particle-particle channels, respec- 
tively. 

In the present studies we consider external fields of the 
pair addition and removal types 



Tz-ext 
pair 



/ drdl-^^{r)YLM{r){P^ir) + P,(r)) 
+ /drS-'^Jr)riM(r)(PJ(r) - P,(r)) (7) 

in addition to the usual external field of the particle-hole 



type 



V-cxt 
ph — 



drJ2v°a'^^L{r)YLM{r)pq{j 

n 



(8) 



Note that the linear response equation ^ formulated in 
the present work contains v^\]^ (r) representing the pair- 
addition/removal fields. 

For the unperturbed response function Rq „]^ appear- 
ing in the linear response equation (JG]) , we adopt the spec- 
tral representation 



Ro^qLir^r',u) 



- E 

nlj.n'l' j' 



{i'j'\\yLm^ 

2L + 1 

{0\A\nlj,n'l'j') (nlj,n'l'j'\B\0) , [a\B\nlj,n' I' f) ^, {nlj.n' I' j'\A\o) ^ 



huj+ie — Enij —E^ 



hie+-B„,j+B„ 



(9) 



r 



using 
HFB 

a(i) 



i't^nlj 



with 

(0\A\nlj,n'l'f)^ ^^l^,^,{r)A<lynij{r), 

{nlj,n'l'j'\A\0)^ ^4>l^{r)A^,,,i'rir), 

the discretized quasiparticle eigenstates of the 
mean-field Hamiltonian. Here 4>nij{i') ~ 

'),(l>niAr))'^ is the two-component radial wave 
function of a quasiparticle state with the quantum num- 
ber nlj, and A and B are 2x2 matrices defined for the 
three kinds of density pq and Pq,±, Eqs.([2]) and ([3]) (see 
Ref.'54|, for details), e is a smearing parameter, for which 
in the present work we adopt e = 0.05 MeV correspond- 
ing to a convolution of the strength function with the 
Lorentzian function of FWHM=100 keV. In the previ- 
ous works ^17i. 48, 49, 54] we used the formalism which 
implements the continuum quasiparticle states satisfy- 
ing the proper boundary conditions connecting to out- 
going waves. The continuum QRPA calculation based 
on this formalism is useful for the states above the par- 
ticle threshold energy, but, on the other hand, it costs 



quite severely a large computation time, especially when 
we choose a small value for the smearing parameter e, 
and it is too demanding for systematic studies. The cal- 
culated lowest 2"*" mode obtained in the present study 
are all located below the neutron separation energy, and 
we use a large box size iimax = 20 fm. In the present 
study, therefore, we adopt the spectrum representation 
using the discretized quasiparticle states obtained in the 
HFB calculation. An argumentation on this point will 
be given later. 



B. Strength function and transition density for 
pair transfer modes 

The HFB-I-QRPA model enables us to calculate the 
transition matrix elements for pair transfer operators. 
Here by the pair transfer operator we mean a generalized 
one-body operator which either add or remove a pair of 
singlet neutrons. We consider the following pair-addition 



operator with multipolarity L, 



- drYLM{f)I{r)P\r) 



and the strength of the pair-removal operator for the in- 
dividual QRPA solution 



- drYLM{f)f{r)^j\rl)ij^{r\), (10) 



LM 



which is essentially the same as the m ultip ole pai r oper- 
ators used in the traditional works [30, IMI, Eo, EH- Here 
and hereafter we omit the isospin index q for simplicity 
since we discuss only the neutron pair transfer and the 
neutron pairing in the present paper. Note that /(r) is a 
form factor of the pair-addition operator, and we choose 
/(r) = 1 in the present analysis. This is a simplifica- 
tion, but even this simple choice is useful for discussions 
on isotopic trends and roles of the surface enhanced pair 
interactions as we show below. 

The strength function for the multipole pair-addition 
operator P^a/ i^ given by 

Sp^ALiE) =Y,5{E - E,l)\{^,lm\pIm\'^o)\\ (11) 

where ^I^o and '^iLM are the ground and the excited states 
obtained in the QRPA, and En, is the QRPA excitation 
energy. In the linear response formalism of QRPA, this 
quantity can be calculated as 

<S'padL(^w) = Im / dr-{Sp+^qL{r,uj)Sp-^qL{r,io)) 

(12) 
using the solutions of Eq. © obtained with the external 
field 



B(PrmL) = J] |(*,l-m|-Plm|*o)I' 

- (2L + l)|(*,io|^io|*o)|'. 
The former is calculated as 



(17) 



SprniL{riUJ) = Im 



dr-{Sp+,qL{r,uj)+Sp-^qL{r,uj)) 



(18) 
under the external field of the multipole pair- removal op- 
erator Plq: 



yext(^) 



(19) 



When we calculate the strength i?(QL) = (2L -|- 
l)|(^iLo|QLo|^o)|^ for the multipole operator of the 
particle-hole type Qlm — J drYLM{f)r^ Pq{r) , we choose 

the external field as V^'^ir) = (r^,0,0)^. 

We also calculate the transition density associated with 

. . pt 

the pair-addition ^o — ^ ^iLM and pair-removal trans- 

P 
fer ^0 — > ^iL-M- The transition density associated 

with the addition/removal of a neutron singlet pair at 

the position r is defined and factorized as 



P^^ir) ^ {^.LM\i'Hri),p\rt)\^o) 



V<^^^r) 



<Lir) \ / 



(13) 



-,qL 



(r) 



which corresponds to an external field identified to the 
pair-addition operator, i.e., V°^^^ = pIq. We evaluate the 
strength 



= {^nM\PHr)\^o)=Y^,Ar)Prir)m 
P^Z]{r) ^ (*.LM|V'(rt)V^(r-i)|*o) 

= {^^LM\P{r)\^,) = Yl,,{r)Pt\r).{2l) 

Here the radial transition densities are expressed using 
the solution of the linear response equation as 



B(PadL) - 5]|(*,lm|p1mI*o)| = 



M 



(2L + l)|(*,io|/'lol*o) 



(14) 



4fM 



C 

27rr2 

C 
27rr2 



Im((5p+^,i,(r,Wj) - 5p-^qL{r,Ui)\^2) 
lm{5p+^qL{r,u:i) + Sp-^qL{r,uJt)\23) 



associated with an individual QRPA solution by integrat- 
ing the strength function S'padi (E) in an energy interval 
E G [Ei — lOe, Ei + lOe], where Ei is the peak energy of 
a peak in SpadhiE), corresponding to a QRPA solution, 
and e is the smearing parameter. An isolated discrete 
peak has the FWHM of 2e = 100 keV. 

Similarly, we can consider the multipole pair-removal 
operator 



where Ei = fiuji is the excitation energy of the QRPA 
solution, and C is a normalization constant. Note that 
the conventional transition density of the particle-hole 
type is given by 



P\ 



(r) = (*,iMl5I^g(^^)^9(^^)l*o) 



Plm ^ j drYl^{f)f{r)^q{r t)^,(r ;) = (Pl^,)t, 



Y, 



ph 



'Mir)P%Lir)^ 



ph 



c 



the associated strength function 



(15) 



pLlW = :^lmSpqLir,uj,) 



(24) 
(25) 



Sp,^,l{E) =J2SiE- -B,L)|(*,:L-M|PLAf |*o)|', (16) 



iM 



The normalization factor C is fixed so that it is consis- 
tent with the strengths B{QL) or B(Pa.dL) of the QRPA 
solution under consideration. 
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FIG. 1: The pairing interaction strength V„{r) for neutrons 
in ^^''Sn, plotted as a function of the radial coordinate r for 
the DDDI-bare', the mix and the volume pairing interactions. 
The position of the half density radius is marked with the 
arrow. 
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FIG. 2: The average neutron pairing gap An in the Sn iso- 
topes. The one calculated with the DDDI-bare' (77 = 0.71) 
is plotted with the solid line. For 100 < vl < 114, results 
obtained with the DDDI-bare' but using 77 = 0.75 is also plot- 
ted with the thin solid line and the diamond symbols. The 
ones calculated with the mix pairing and the volume pair- 
ing are plotted with the dotted and the dashed lines, respec- 
tively. The experimental neutron gap estimated u sing the 
three-point formula of the odd-even mass difference[6^ and 
the mass compilation [6q] is also plotted with the circles. 



III. EFFECTIVE PAIRING INTERACTIONS 

In the present paper we intend to disclose influences of 
the surface enhanced pairing on the low-lying quadrupole 
vibrational mode. For this purpose, we compare three 
versions of DDDI, Eq.®, which have different proper- 
ties concerning the surface enhancement. They are 1) 
the density-dependent pairing interaction which has a 
connection to the bare nuclear force via the scattering 
length, 2) the density-independent delta interaction, and 
3) the so-called mix pairing interaction |60l - [63 | . 

1) The first is the one which we adopt mainly in the 
following analysis, and it is a version of DDDI which 



reproduces properties of the bare nuclear force relevant 
for the pairing, in particular, the scattering length of the 
bare forced [23, [2I [13, [11]. Among such DDDFs, we 
adopt the parametrization 



V„{r) = t;o 1 -77 



Pn(r)\" 



Pc 



(26) 



for the interaction of neutrons [23l |49|. Here the pa- 
rameter vq is chosen to reproduce the scattering length 
a = — 18.5fm in the ^S channel. The density-dependence 
characterized by the power a and and the prefactor t] is 
adjusted to reproduce the neutron pairing gap in dilute 
and uniform matter, more specifically the pairing gap 
obtained with the bare force and the BCS approxima- 
tion. We find a parameter set vq — —458.4 MeV fm~^, 
Pc = 0.08fm~3, a = 0.845, and 77 = 0.845 for the quasi- 
particle energy cut-off with Ecut — 60MeV[2j|. However 
this parameter set leads to neutron pairing gaps in Sn iso- 
topes which are systematically smaller than the empirical 
values estimated from the odd-even mass difference us- 
ing the three-point formula [65J. We therefore vary rj and 
fix to 77 = 0.71 so as to reproduce the empirical gap in 
^■^°Sn[4^] while keeping the scattering length relevant to 
the low-density limit. We call this parameter set DDDI- 
bare' since it is refereed to the bare nucleon force, but 
with a small modification. 

2) The second is the simplest density-independent 
delta interaction; 



V„{r) =vo. 



(27) 



The only relevant parameter vq is fixed to vq ~ —195 
MeVfm^ to reproduce the gap in ^^°Sn, as is done in 
Ref.[60, [63 • The density- independent delta interaction 
is often called the volume pairing inter action [36|- 

3) The third is the so-called mix pairing interaction [6ll 
|62| . which has a linear dependence on the total nucleon 
density p(r), and given by 



Vnir)=vo[l- 



1 fpir] 



Pa 



(28) 



with po = 0.16fm^'^. Here the parameter vq is fixed to re- 
produce the gap in ^^"Sn, leading to vq = —292 MeVfm'^. 
The pairing force strength Vn{r) is compared among 
the three kinds of interaction in FigHJ which plots Vn{r) 
as a function of the radial coordinate in ^'^^Sn. It is 
seen that the DDDI-bare' exhibits the strongest radial 
dependence, i.e., the strongest density dependence, and 
that the interaction strength in the region outer than the 
nuclear surface ( corresponding to the low density limit 
p{r) — i> ) is more than twice larger than that of the vol- 
ume pairing. The mix pairing exhibits an intermediate 
radial- (density-) dependence. Since these three versions 
of DDDI have very different force strength in the external 
region we expect to see the difference in physical observ- 
ables related to the surface region and/or in the exterior 
of nuclei. 
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FIG. 3: (a) The excitation energy of the calculated 2^ QRPA solutions obtained with the DDDI-bare', the mix and the volume 
pairing interactions, plotted with the diamonds connected by the solid line, the dotted line, and the dashed line, respectively. 
The experimental 2^ energv[72| is plotted with the circles, (b) The same but for _B(E2;0ji — >■ 2^). The experimental data 
is taken from Refs. [73 - [77l ] (c) The same but for the neutron quadrupole strength _B(N2;0j; -^ 2^), and (d) for the ratio 
{Mn/Mp)/{N/Z) of the neutron and proton amplitudes associated with the 2+ QRPA solution. 



We compare in Figl2] the calculated average neutron 
pairing gap A„ = J A{r)nPn{r)dr/ J pn{r)dr obtained 
with the three interactions. It is seen that they give 
the average pairing gaps which are almost identical for 
A - 116 - 130. In neutron-rich nuclei A = 134 - 138 
beyond the A^ = 82 shell closure, a small difference of 
the order of < 150 keV is seen (except the subshell nu- 
cleus ^^°Sn where the volume pairing gives the vanishing 
gap). It is pointed out in the literature [OtMijl that the 
surface enhancement of the pairing interaction influences 
the average paring gap only to a small extent unless one 
sees the trend in a very wide range of mass|6lL[6a|. and 
the difference may grow only in very neutron-rich nuclei 
close to the drip-line 60, 63] . On the other hand, it is seen 
that the DDDI-bare' gives an overestimate of the average 
pairing gap for the isotopes with A ^ 116, by about <; 200 
keV compared with the other two interactions. For a fair 
comparison in the following analysis, we modify slightly 
the value of ry, i.e., t] = 0.75, only for A < 116 so that it 
gives the average pairing gap comparable with those of 
the volume and the mixed pairing (cf. the thin solid line 
in Figlll). 



IV. RESULTS 

In the first two subsections we will mainly discuss re- 
sults which we obtain with the effective pairing interac- 
tion DDDI-bare' having the strong increase of the inter- 
action strength in the surface and the external regions. A 
comparative analysis using different pairing interactions 
will be presented in the subsection IIVCI 



A. The lowest-lying quadrupole state 

Let us first discuss the basic properties of the calcu- 
lated lowest-lying quadrupole state, on which we focus 
in the present paper. Figure [3lja) and (b) shows the ex- 
citation energy £^(2^) and B(E2;0+ -^ 2^) of the 2^ 
QRPA solution in the Sn isotopes. The results are simi- 
lar to those in the full self-consistent QRPA calculation 
[68| using the same Skyrme interaction SLy4. It is seen in 
Fig. [3);a) and (b) that the calculated E{2+) and B{E2) 
reproduce reasonably well the experimental values. In 
the present QRPA calculation the lowest quadrupole 
mode shows an instability in the A = 110 — 114 and 
A > 144 isotopes as is also seen in Ref . [68| . We do not 
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FIG. 4: The transition densities associated with the 2^ QRPA 
solution (a) in ^^''Sn and (b) in ""Sn. The soUd and the dot- 
ted curves are the particle-hole transition density r^p*'''(r) 
for neutrons and protons respectively while the dashed and 
the dotdashed curves represent the pair-addition transition 
density r^P^'^'^'{r) and the pair-removal transition density 
^2p{rm)|-^-j of neutrons, respectively. The arrow indicates the 



half density radius Ri/2- 
the DDDI-bare'. 



The effective pairing interaction is 



discuss these cases and the A = 116, 144 cases where 
the QRPA excitation energy is less than 0.5 MeV since 
the quadrupole mode then needs to be described as an 
anharmonic vibration, which is beyond the scope of the 
present work. Note that there exist a few other QRPA 
calculations for the neutron-rich Sn[69ll70|. 

The lowest-lying quadrupole state in isotopes beyond 
the iV = 82 shell closure are of special interest. Figure 
^c) shows the neutron quadrupole strength B(N2; 0+ -> 
2^) calculated in the same way as i?(E2), but for neu- 
trons instead of protons. The isotopic dependences of 
S(E2) and B(N2) look similar for N < 82, but the values 
of B(N2) are apparently larger than the proton strength 
B(E2)/e^ as the neutron number exceeds the magic num- 
ber N — 82. This behavior is s een clearly in the ratio 
Mn/Mp = ^B(N2)/(B(E2)/e2) of the neutron and pro- 
ton amplitudes, plotted in Fig. ^d). If the low- lying 
quadrupole mode is a collective surface oscillation of liq- 
uid drop of uniform mixture of neutrons and protons, the 
ratio Mn/Mp would be N/Z. The calculated double ra- 



tio {Mn/Mp)/{N/Z) is around 1.2-1.3 for A < 132, being 
not very different from the liquid drop limit 1, and it is 
consistent with the typical values (« 1.3 — 1.5) evaluated 
experimentally in stable isotopes ^^^^^^^Sn[78|. The dou- 
ble ratio, however, increases suddenly in A > 132, and 
amounts to 1.6-1.7. The ratio Mn/Mp itself reaches 2.7- 
3.1. In the isotopes beyond N = 82, the amplitude Af„ 
of neutrons dominate over that Mp of protons in the first 
2+ state. Note that the {Mn/Mp)/{N/Z) ratio obtained 
in Ref.jZO, UM i^ somewhat smaller than our values by 
about 0.2. 

A characteristic feature of the lowest-lying quadrupole 
states beyond the iV = 82 shell gap can be demon- 
strated more explicitly in terms of the transition den- 
sities, which are shown in Fig|3] for (a) ^■^^Sn and (b) 
^^"Sn. In ^^'^Sn, the transition densities exhibit the fea- 
tures typical of the low-lying quadrupole surface vibra- 
tion. Namely, the particle-hole transition densities for 
neutrons and protons, and the pair-addition and -removal 
transition densities of neutrons exhibit similar radial pro- 
file all of which have peaks around the nuclear surface 
r ^ Ri/2 ~ 5.5 fm {R1/2 denotes the half-density ra- 
dius defined by p{r) — 0.08 fm^'^). In ^'^^Sn, in contrast, 
the peak positions of the neutron and proton particle- 
hole transition density Pn\{r) deviates significantly with 
each other: the neutron transition density has amplitude 
which extends toward the exterior of the nucleus while 
the proton transition density does not. 

The transition density p(^'^)(r) of the neutron pair- 
addition mode is most characteristic in ^'^'*Sn. It has the 
largest amplitude in the external region r^Ri/2 (~ 5.7 
fm) among the three kinds of transition densities. The 
position where the amplitude of pt^'^) (r) is maximum, 
r w 6.6 fm, deviates by about Ifm from the nuclear sur- 
face (r « 5.7 fm), and it is even outside of the peak po- 
sition of the neutron particle-hole transition density. In 
addition, it is extended significantly toward the outside 
of the nucleus, and the tail reaches far to 12-13 fm (in the 
scale of the figure) . The extension of the amplitude in the 
external region is much more significant than that in the 
neutron particle-hole transition density fP^{r). We note 
also that the transition density p(''™)(r) of the neutron 
pair-removal mode is very small in comparison with the 
pair-addition transition density P'^^'^\r)^ and the domi- 
nance of the pair-addition amplitude is commonly seen 
in aU the cases of ^'^^"^'^^Sn. 



B. Quadrupole pair transfer 

Figure [5] shows the strength functions S'pad2 {E) and 
'S'piin2 (£-) for the quadrupole pair-addition and -removal 
modes in ^■^''Sn and ^^"Sn. The pair-addition strength 
associated with the first 2+ mode in ^'^'*Sn is signifi- 
cantly larger than that in -"^^"Sn, and it reflects the large 
pair-addition amplitude seen in FigH) The unperturbed 
strength functions are also plotted for comparison (the 
dotted curve). The first peak in the unperturbed cal- 




3 4 

E (MeV) 

FIG. 5: (a) The transition strength function Spad2{E) for the quadrupole pair-addition mode in Sn, plotted with the solid 
curve. The thin dashed curve is the result of the continuum QRPA calculation while the dotted curve represents the unperturbed 
response for the two-quasiparticle excitations. The effective pairing interaction is the DDDI-bare'. The smearing parameter is 
e = 50 keV. The arrows indicate the calculated one- and two-neutron separation energies Si„ and 82,1- (b) The same as (a) 
but for the transition strength function Sprni2{E) for the quadrupole pair-removal mode. (c)(d) The same as (a) and (b), but 
for ^20 Sn. 



culation of "'^^^Sn corresponds to the transition to the 
two-quasi-neutron configuration [l/7/2]/=2- The pair- 
addition strength for the 2^ state in the full calculation 
is about ten times larger than that in the unperturbed 
calculation, thus it indicates a significant enhancement 
due to a collectivity (the RPA correlation in the other 
words) in the pair-addition mode. Another feature ob- 
served in ^'^'*Sn is that the RPA correlation effect is large 
not only for the lowest mode but also for higher lying 
modes: the second peak of Spad2{E) also has a large 
pair-addition strength comparable with the first peak, 
and there exist several peaks in the excitation energy re- 
gion E^A MeV with significant pair-addition strength. 
The peaks above the one-neutron separation energy Sin 
are discretized continuum states. The continuum QRPA 
calculation [4^ [5J| is appropriate for the excitation en- 
ergies above Sin- In fact, the result of the continuum 
QRPA calculation, shown with the thin dashed curve in 
Fig. [SJ indicates that the proper treatment of the contin- 
uum modifies the strength function for E > Sin- On the 



other hand, it has only a minor effect on the strengths 
below Sin, and it does not affect the discussions below. 
Note also that the pair-removal strength in ^^''^Sn is ex- 
tremely small in the displayed energy range. 



1. pair- addition mode 

We shall discuss properties of the pair-addition mode 
in detail. Figure HJa) is the systematics of the pair- 
addition strength i3(Pad2) associated with the 2^ mode. 
A feature immediately noticed in the figure is a strong 
variation along the isotopic chain. i3(Pad2) takes a 
rather small value for A < 132, but it shows a big jump 
a,t A — 132, where the strength is larger than those in 
any other isotopes with A < 132. Proceeding from ^"^^Sn 
to^-^^Sn, ^'^^Sn and ^'^^Sn, the strength decreases gradu- 
ally while keeping relatively large values. It then turns to 
increase in ^^'^Sn. and takes a very large value in ^^^Sn. 

Figure[7Ja) shows the pair-addition transition densities 




100 



no 



120 . 130 
A 



140 



FIG. 6: (a) The calculated transition strength _B(Pad2) of 
the neutron pair-addition transfer feeding the 2f states in Sn 
isotopes. The diamonds connected with the solid line are the 
results obtained with DDDI-bare' while the dotted and the 
dashed lines (the circles and the crosses) are for the mix and 
the volume pairing interactions, respectively, (b) The same 
but for the strength _B(Prm2) of the quadrupole pair removal 
transfer. 



p(ad)^y.^ associated with the lowest quadrupole mode in 
the isotopes ^^^"^''^Sn beyond the A^ = 82 shell closure. 
It is seen that they all exhibit sizable amplitudes extend- 
ing to r ^ 8 fm as is found in Fig|4]in the case of ^"^^Sn. 
The spatial extension of the pair-addition transition den- 
sity develops further in ^*°Sn and ^^^Sn, for which we see 
significant amplitude even around r ^ 12—13 fm, which is 
about '^ 7 f m apart from the half density surface R1/2 = 
5.6 — 5.8 fm. This is contrasted with the pair-addition 
transition density for the isotopes i20-i30gj^ below the 
shell closure (displayed in Figl^b)). In these nuclei spa- 
tial extension is hardly seen in the pair-addition transi- 
tion density, and the transition densities is mostly con- 



centrated in the internal and surface region r <i Ri 



/2 



a 8 fm. The spatial extension emerges suddenly at 
132. Comparing the pair-addition transition density 



fm 
A = 

in ^^^Sn (the largest one in Figlll^b)) with that in ^^^Sn 
(the third largest one among the solid curves in FiglTl^a)), 
the maximum amplitude is almost the same, but the pair- 
addition strength S(Pad2) in "6Sn (S(Pad2) = 1.42) is 
about twice as large as that in ^^^Sn (i3(Pad2) = 0.82) 
because of the spatial extension in ^'^^Sn. 



The sudden change across the N = 82 shell closure may 
be ascribed to the fact that the isotopes with A > 132 
have a feature of weakly bound nuclei although they 
are located far from the neutron drip-line. The cal- 
culated one-neutron separation energy suddenly change 
from 5i„ = 7.62 MeV in ^^sgn to Sin = 3.32 MeV in 
^■^^Sn. This arises from the shell gap at iV = 82 and 
from the fact that the neutron single-particle orbits I/7/2 
and 3p3/2 lying above the TV = 82 shell gap are bound 
only weakly ( the Hartree-Fock single-particle energies 



in ■^ Sn are ei/^ 



/2 



-1.99 MeV and 63^. 



'3/2 



-0.25 

MeV). These orbits have significantly small binding en- 
ergies compared with the orbit /111/2 just below the shell 
gap (e;ijj^2 — —7.68 MeV). Since the quadrupole pair- 
addition mode in ^^^Sn and heavier isotopes necessarily 
contain two-quasiparticle components consisting of these 
weakly bound orbits, the pair-addition transition density 
should inherit the long tail arising from the weak bind- 
ing. In particular, approaching ^"^^Sn and ^^^Sn, neu- 
trons start to occupy the 3p3/2 orbit which has a very 
small binding energy. 

We remark that the closing of shells has a strong in- 
fluence on the pair transfer in general 13 11]. For instance, 
the neutron shell is closed at ^ = 132 {N = 82) and 
the calculated neutron pairing gap vanishes there. In 
this case, the QRPA equation ([6]) is decomposed into 
the particle-hole, the particle-particle, and the hole-hole 
sectors. The pair-addition mode in ^'^^Sn belongs to 
the particle-particle sector, and the lowest-lying QRPA 
solution in this sector has a character of a correlated 
pair of neutrons built upon the ^'^^Sn core. It is noth- 
ing but the quadrupole pair vibrational mode in closed 
shell nuclei, and in this case a large collectivity in the 
pair-addition mode is known, with a typical example of 
^°^Pb[3li|40,|il|. This is, however, only one of the mech- 
anisms of the large pair-addition strength i3(Pad2) seen 
at ^ = 132 (Af = 82). We note that the strength in "2Sn 
for the DDDI-bare' interaction is more than twice larger 
than that in -'^™Sn (another shell closure at iV = 50), 
and this relative enhancement at ^'^^Sn is a consequence 
of the characteristic spatial extension of the transition 
density in the isotopes beyond the N = 82 shell closure. 



2. pair-removal mode 

Figure [6] (b) shows the systematics of the pair-removal 
strength i3(Prm2) for transition from the ground state 
of an isotope A to the 2^ state in the isotope A — 2. It 
is seen that the pair-removal strength is generally small, 
compared with the pair-addition strength. Another trend 
seen here is that the pair-removal strength is very small 
at the beginning of shells (at A sa 102, 134) and of a 
subshell (at A — 142); it takes a value comparable with 
the pair-addition strength B(Pad2) around the mid-shell 
{A ^ 120), and it is relatively large at ^ = 100, 132 
because of the shell closures at iV = 50 and A^ = 82 in 
parallel to the situation in •^°^Pb.j31[ 
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FIG. 7: (a) The transition density r^P^'^'^'{r) for the neutron pair-addition transfer associated with the 2f QRPA solutions, 



-132 



'Sn. The pair-addition transition strength B(Pad2) is used to normahze, and the phase is 

Sn while the dotted curve is used 

""Sn. (c) The transition 

^Sn. The 



calculated in the isotopes 

chosen so that the transition matrix element becomes positive. The solid curve is for ^ 

for ^*''~^'*^Sn. The effective pairing interaction is the DDDI-bare'. (b) The same as (a), but for 

density r'^ P^'^'^' (r) for the neutron pair-removal transfer associated with the 2^ QRPA solutions, calculated in 

solid and the dotted curves are for ^^*^^'^^Sn and ^^"'^^'^Sn, respectively. 



Figure W[c) displays the pair- removal transition densi- 
ties p(''™) (r) associated with the 2'^ mode in the isotopes 
from ^-^^Sn (located in the middle of the shell) to ^'^■^Sn 
(at the end of the shell). It is seen that the transition 
densities p(''°) (r) of the pair-removal mode in these iso- 
topes does not exhibit the spatial extension toward the 
outside, rather it is concentrated only around and inside 
the surface region r ^ R1/2 + 2 f m ~ 8 fm. This feature 
originates from the fact that the pair-removal transition 
consists mainly of the obits situated near the Fermi en- 
ergy and those bound more deeply. In the isotopes with 
A < 132, these are strongly bound orbits with the bind- 
ing energy e^.p. ^ -7 MeV (e.g, Ch^^^^ = -7.68 MeV in 
"'^^^Sn) and they do not have tails extending far outside 
the nuclear surface. 



C. Sensitivity to the surface-enhanced pairing 

We shall now discuss the role of the residual pairing 
interaction and the sensitivity of the quadrupole pair 
transfer to the enhancement of the pairing interaction 



in the surface and the external regions of a nucleus. In 
FiglHlwe compare the pair-addition strength function in 
^^''Sn and ^^°Sn, calculated with two pairing interactions, 
the DDDI-bare' and the volume pairing, which have very 
different density-dependence in the pairing interaction 
strength (cf. FigH]). In the case of ^^''Sn we see a 
large difference in the quadrupole pair-addition strength 
in ^'^^Sn while in the case of ^^°Sn the difference is small. 

A systematic comparison is made in Fig [5] for the pair- 
addition and -removal strengths i3(Pad2) and i?(Prm2) 
of the 2^ mode calculated with the three kinds of the 
pairing interactions, including the mix pairing as well. It 
is seen that the pair-addition strength exhibits a sensi- 
tivity in all the nuclei to some extent, and the sensitivity 
becomes prominent in ^'^^~^'*^Sn. The ratio of i3(Pad2) 
obtained with the DDDI-bare' and that with the volume 
pairing is 2.3, 2.8, 4.0, 3.5, 2.8 for A = 134,136,138, 140, 
142, respectively, while for A < 132 it is in the range 
1.4-2.0. Note, on the other hand, that the pair-removal 
strength i?(Prm2) is much less sensitive in all the iso- 
topes. The difference in B(Prm2) for the different pairing 
interactions is negligible in many isotopes. 
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FIG. 8: Comparison of the transition strength functions 
<S'pad2(S) for the neutron quadrupole pair-addition transfer, 
calculated with the volume pairing interaction (solid curve) 
and with the DDDl-bare' (dotted curve), (a) and (b) are for 
^^"Sn and ^^°Sn, respectively. 



Let us analyze the mechanism of the sensitivity of the 
pair-addition strength by looking into the transition den- 
sities. We shovif in Fig. [9] (a)-(c) the transition density 
P^^^^ir) of the pair-addition mode calculated with dif- 
ferent pairing interactions in the three representative Sn 
isotopes. In this figure we see an intimate relation be- 
tween the degree of sensitivity and the radial profile of the 
transition density. Namely a large sensitivity observed in 
^■^^Sn and ^^^Sn can be linked to the pair- addition tran- 
sition density extending far outside the nuclear surface, 
r^Ri/2 + 2fm. The sensitivity is small in ^^"Sn where 
the transition density is confined mostly within the sur- 
face region r ^ R1/2 + 2fm. 

For more detailed discussion let us examine the values 
of the pair-addition transition density p'^^){r) in ^'^^Sn 
at r = 5.8,7.8,9.8 fm, corresponding approximately to 
the half density radius r ~ R1/2 (= 5.70 fm), and r — 
Ri/2 + 2 fm, and r = R1/2 + 4 fm. The ratio of p('"^'>{r) 
obtained with the DDDI-bare' and that with the volume 
pairing is 1.45, 1.82 and 2.71 for r = 5.8,7.8,9.8 fm, 
respectively, indicating an increase of the sensitivity as 



moving from the surface toward the outside. 

It is also useful to compare p('"^)(r) with the tran- 
sition density Punp (r) associated with the unperturbed 
pair-addition transfer (i.e. the transition associated with 
the lowest-energy two-quasiparticle excitation obtained 
by neglecting the residual pairing interaction), which is 
also shown in Fig |9] with the dotted curve. The ratio be- 
tween the transition density p('*"^)(r) obtained with the 
DDDI-bare' (or the volume pairing) and the unperturbed 

transition density Punp (r) is a measure of the infiuence 
of the RPA correlation (the configuration mixing) caused 
by the residual pairing interaction. In the case of the 

DDDI-bare', the ratio P('''^'('')/^unp^(^) reads 2.7 at the 
nuclear surface r = R1/2, and it increases to 3.8-6.2 at 
r = {R^/2 + 2fm) - (P1/2 + 4fm) w 8 - lOfm. This is in 
contrast to the volume pairing causing much smaller ratio 
1.9-2.1-2.3 at the same positions. We see even stronger 
sensitivity to the density-dependence of the pair interac- 
tion in ^''^Sn, for which the tail extends further outside. 

We thus see that the strong sensitivity of the pair- 
addition transfer originates from the RPA correlation 
taking place in the surface and exterior regions in ^^"^Sn 
and ^''^Sn. We recall the reader here that the three 
pairing interactions have interaction strengths different 
by a factor of more than 2 (cf. FigH]) at the positions 
r <:, Pi /2 + 2 fm corresponding to low densities p/po ^ 0.1. 

A contrasting case is ^^°Sn shown in Fig.[5fc), in which 
we see very little the RPA correlation effect and the dif- 
ference between the DDDI-bare' and the volume pairing 
interactions. This is what we can expect from the dis- 
cussions above as there is essentially no significant spatial 
extension of the pair-addition transition density in this 
tightly bound stable isotope, as seen in Fig.[9jc). 

In Fig. m^d) we show the transition density p(''™)(r) 
of the pair- removal mode in ^^^Sn. In this case the pair- 
removal amplitude is concentrated only in the vicinity 
and the inside of the surface r ^ R1/2 + 2 fm, exhibiting a 
very different radial profile from that of the pair-addition 
mode in "4,i42gjj ^pig_ [gj^^) ^^^ ^y^y j^ -g ^jg^ ^^^^ 

that the pair-removal transition density p(''™)(r) does 
not depend on which we adopt the DDDI-bare' or the 
volume pairing. We can interpret this insensitivity in 
terms of Fig. [TJ the effective pairing strengths of the 
three interactions are comparable in size on average in 
the region \r — R1/2 1 ^ 2 f m where the transition strength 
is concentrated. 

Finally we remark that the three pairing interactions 
give essentially identical results for i?(2^) and P(E2) 
except for A > 140, as is seen in Fig. [3] (a) and (b). 
The excitation energy and the P(E2) of the 2f state 
in superfiuid nuclei is known to depend on the pairing 
correlation [l|. However, the results indicate that these 
quantities are not influenced very much by the density- 
dependence of the pairing interaction as far as the aver- 
age pairing gap is constrained to the same value. 
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FIG. 9: Effects of the pairing residual interactions on the neutron pair-addition transition density r''P^'^^'{r) feeding the 2\ 
states, calculated for the isotopes (a) ^*^Sn, (b) ^^*Sn, and (c) ^^"Sn. The solid and the dashed curves are the results obtained 
with the DDDI-bare' and the volume pairing, respectively, while the unperturbed pair-addition transition density r^Punp (r) 
neglecting the residual interaction is also plotted with the dotted curve. The arrow indicates the position of the half density 
radius i?i/2. (d)The same as (a-c), but for the pair-removal transition density r^P^'^^'ir) of the 2\ pair-removal mode in ^^^Sn. 



V. CONCLUSIONS 

Wc have investigated the neutron pair transfer modes 
with the quadrupole multipolarity feeding the low-lying 
quadrupole states in the neutron-rich Sn isotopes by 
means of the QRPA based on the Skyrme-Hartree-Fock- 
Bogoliubov mean-fields. As a model of the surface- 
enhanced paring, we used the density-dependent delta 
interaction (DDDI) whose strength increases at low nu- 
cleon densities reflecting the large scattering length of the 
^S bare interaction in the free space. 

The numerical calculations indicate that the lowest- 
lying quadrupole state exhibits a sudden change of its 
character as the neutron number increases across the 
magic number N = 82, beyond which features of weakly 
binding emerge in the neutron motion. The most notice- 
able feature is found in the pair-addition mode. Namely 
the transition density of the pair-addition mode feed- 
ing the 2^ state extends spatially far outside the nu- 
clear surface with a long tail reaching r ~ 12 — 13 fm 
in the isotopes ^^'^^Sn. Accordingly the strength of the 
quadrupole neutron pair-addition transition is enhanced 
for A > 132. The enhanced pair-addition transfer is a 



consequence of the RPA correlation caused by the pair- 
ing interaction acting strongly in the surface and external 
regions of a nucleus. 

Comparing with other DDDFs having different density 
dependence, i.e. the density-independent volume pairing 
and the mixed pairing having an intermediate density de- 
pendence, it is found that the quadrupole pair-addition 
mode exhibits generally sensitivity to the increase of the 
pairing interaction strength in the surface and the ex- 
ternal regions of a nucleus. The sensitivity of the pair- 
addition mode to the three interactions is about a factor 
of 1.4-2.0 in the isotopes with A < 132, and it is mag- 
nified to a factor of 2.0-4.0 in the isotopes A > 132. 
Therefore one can expect that the neutron pair-addition 
transfer feeding the 2'^ state in the neighboring A + 2 
isotope, especially those feeding the isotopes A > 132, 
provides a probe to the enhancement of the pairing in- 
teraction strength expected in the surface and the ex- 
ternal regions. On the other hand, the neutron pair- 
removal transfer feeding the 2^ state in the A — 2 isotope 
does not have the sensitivity, and neither B{E2) nor the 
quadrupole neutron strength B{N2) does, as far as the 
Sn isotopes analyzed in the present work are concerned. 



13 



Acknowledgments 

One of the authors (M.M.) thanks Prof. M. Kawai, 
Prof. M. Ichimura and Prof. A. Vitturi for useful 
discussions. A part of the numerical calculations were 
performed on the NEC SX-8 supercomputer systems at 
Yukawa Institute for Theoretical Physics, Kyoto Univer- 



sity, and at Research Center for Nuclear Physics, Osaka 
University. The work was supported by the Grant-in- 
Aid for Scientific Research (Nos. 20540259, 21105507 
and 21340073) from the Japan Society for the Promo- 
tion of Science, and also by the JSPS Core-to-Core Pro- 
gram, International Research Network for Exotic Femto 
Systems(EFES). 



[1] A. Bohr and B. R. Mottelson, Nuclear Structure vol. II [27 

(Benjamin, 1975). 

[2] D. M. Brink and R. A. Broglia, Nuclear Superfluidity: [28 

Pairing in Finite Systems (Cambridge University Press, 

Cambridge, 2005). [29 

[3] Y. R. Shimizu, J. D. Garrett, R. A. Broglia, M. Gallardo, [30 

and E. Vigezzi, Rev. Mod. Phys. 61, 131 (1989). [31 

[4] D. J. Dean and M. Hjorth- Jensen, Rev. Mod. Phys. 75, 

607 (2003). 
[5] G. F. Bertsch and H. Esbensen, Ann. Phys. 209, 327 

(1991); H. Esbensen and G. F. Bertsch, Nucl. Phys. [32 

A542, 310 (1992). [33 

[6] K. Ikeda, Nucl. Phys. A538, 355c (1992). [34 

[7] M. V. Zhukov, B. V. Danilin, D. V. Fedorov, J. M. Bang, [35 

I. J. Thompson, and J. S. Vaagen, Phys. Rep. 231, 151 [36 

(1993). 

[8] K. Hagino and H. Sagawa, Phys. Rev. C 72, 044321 

(2005). [37 

[9] K. Hagino, H. Sagawa, J. Carbonell, and P. Schuck, Phys. 

Rev. Lett. 99, 022506 (2007). [38 

[10] Y. Kikuchi, K. Kato, T. Myo, M. Takashina, and 

K. Ikeda, Phys. Rev. C 81, 044308 (2010). [39 

[11] T. Nakamura et al., Phys. Rev. Lett. 96, 252502 (2006). 
[12] P. Mueller et al., Phys. Rev. Lett. 99, 252501 (2007). [40 

[13] R. H. Ibarra, N. Austern, M. Vallieres, and D. H. Feng, [41 

Nucl. Phys. A288, 397 (1977). 

[14] F. A. Janouch and R. J. Liotta, Phys. Rev. C 27, 896 [42 

(1983). [43 

[15] F. Catara, A. Insolia, E. Maglione, and A. Vitturi, Phys. [44 

Rev. C 29, 1091 (1984). [45 

[16] L. Ferreira, R. Liotta, C. H. Dasso, R. A. Broglia, and [46 

A. Winther, Nucl. Phys. A426, 276 (1984). [47 

[17] M. Matsuo, K. Mizuyama, and Y. Serizawa, Phys. Rev. [48 

C 71, 064326 (2005). 
[18] N. Pillet, N. Sandulescu, and P. Schuck, Phys. Rev. C [49 

76, 024310 (2007). 
[19] N. Pillet, N. Sandulescu, P. Schuck, and J. -F. Berger, [50 

Phys. Rev. C 81, 034307 (2010). 
[20] A. J. Leggett, in Modern Trends in the Theory of Con- [51 

densed Matter (Lecture Note in Physics vol 115), edited 

by A. Pekalski and R. Przystawa, (Berlin Springer- [52 

Verlag, 1980). 

[21] A. J. Leggett, J. de Phys. 41, C7-19 (1980). [53 

[22] P. Nozieres and S. Schmitt-Rink, J. Low Temp. Phys. 59, 

195 (1985). [54 

[23] M. Matsuo, Phys. Rev. C 73, 044309 (2006). 
[24] J. Margueron, H. Sagawa, and K. Hagino, Phys. Rev. C [55 

77, 054309 (2008). 
[25] T. Abe and R. Seki, Phys. Rev. C 79, 054002 (2009). [56 
[26] S. S. Pankratov, E. E. Saperstein, M. V. Zverev, 

M. Baldo, and U. Lombardo, Phys. Rev. C 79, 024309 
(2009). [57 



Y. Kanada-En'yo, N. Hinohara, T. Suhara, and 

P. Schuck, Phys. Rev. C 79, 054305 (2009). 

A. Pastore, F. Barranco, R. A. Broglia and E. Vigezzi, 

Phys. Rev. C 78, 024315 (2008). 

S. Yoshida, Nucl. Phys. 33, 685 (1962). 

D. R. Bes and R. A. BrogUa, Nucl. Phys. 80, 289 (1966). 
R. A. Broglia, O. Hansen, and C. Riedel, "Two-Neutron 
Transfer Reactions and the Pairing Model" in Advances 
in Nuclear Physics vol.6, (Plenum, New York, 1973) ed. 
by M. Baranger and E. Vogt, pp.287-457. 

A. Chatterjee et al, Phys. Rev. Lett. 101, 032701 (2008). 
I. Tanihata et al, Phys. Rev. Lett. 100, 192502 (2008). 
N. Keeley et al, Phys. Lett. B646, 222 (2007). 

M. S. Golovkov et al, Phys. Lett. B672, 22 (2009). 
J. Dobaczewski, W. Nazarewicz, T. R. Werner, 
J. F. Berger, C. R. Chinn, and J. Decharge, Phys. Rev. 
C 53, 2809 (1996). 

E. Khan, N. Sadulescu, N. Van Giai, and M. Grasso, 
Phys. Rev. C 69, 014314 (2004). 

B. Avez, C. Simenel, and Ph. Chomaz, Phys. Rev. C 78, 
044318 (2008). 

E. Khan, M. Grasso, and J. Margueron, Phys. Rev. C 
80, 044328 (2009). 

D. R. Bes and R. A. Broglia, Phys. Rev. C 3, 2349 (1971). 
R. A. Brogha, D. R. Bes, and B. S. Nilsson, Phys. Lett. 
B50, 213 (1974). 

H. J. Ong et al, Phys. Rev. C 73, 024610 (2006). 
H. J. Ong et al, Phys. Rev. C 78, 014308 (2008). 
Z. Elekes et al, Phys. Rev. C 79, 011302 (2009). 
J. K. Jewell et al, Phys. Lett. B454, 181 (1999). 

E. Khan et al, Phys. Lett. B490, 45 (2000). 

N. Iwasa et al, Phys. Rev. C 78, 024306 (2008). 

Y. Serizawa and M. Matsuo, Prog. Theor. Phys. 121, 97 

(2009). 

M. Matsuo, Y. Serizawa, and K. Mizuyama, Nucl. Phys. 

A788, 307c (2007). 

N. Paar, D. Vretenar, E. Khan, and G. Colo, 

Rep. Prog. Phys. 70, 691 (2007). 

M. Igarashi, K.-I. Kubo, and K. Yagi, Phys. Rep. 199, 1 

(1991). 

W. von Oertzen and A. Vitturi, Rep. Prog. Phys. 64, 

1247 (2001). 

G. Potel, B. F. Bayman, F. Brancco, E. Vigezzi, and 

R. A. Brogha, arXiv:0906.4298 

M. Matsuo, Nucl. Phys. A 696, 371 (2001); Prog. Theor. 

Phys. Suppl. 146, 110 (2002). 

M. Bender, P.-H. Heenen and P.-G. Reinhard, Rev. Mod. 

Phys. 75, 121 (2003). 

E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and 

R. Schaeffer, Nucl. Phys. A635, 231 (1998); Nucl. Phys. 

A 643, 441 (1998). 

E. Garrido, P. Sarriguren, E. Moya de Guerra, and 



14 



P. Schuck, Phys. Rev. C 60, 064312 (1999). 
[58] E. Garrido, P. Sarriguren, E. Moya de Guerra, U. Lom- 

bardo, P. Schuck, and H. J. Schulze, Phys. Rev. C 63, 

037304 (2001). 
[59] R. R. Chasman, Phys. Rev. C 14, 1935 (1976). 
[60] J. Dobaczewski, W. Nazarewicz, and P.-G. Reinhard, 

NucL Phys. A693, 361 (2001). 
[61] J. Dobaczewski and W. Nazarewicz, Prog. Theor. Phys. 

Suppl. 146, 70 (2002). 
[62] J. Dobaczewski, W. Nazarewicz, and M. V. Stoitsov, 

Euro. Phys. J. A15, 21 (2002). 
[63] M. Yamagami, Phys. Rev. C 72, 064308 (2005). 
[64] M. Yamagami, Y. R. Shimizu, and T. Nakatsukasa, Phys. 

Rev. C 80, 064301 (2009). 
[65] W. Satula, J. Dobaczewski, and W. Nazarewicz, Phys. 

Rev. Lett. 81, 3599 (1998). 
[66] A. H. Wapstra, G. Audi, and C. Thibault, Nucl. Phys. 

A729, 129 (2003). 



[67] J. Terasaki, H. Flocard, P.-H. Heenen, and P. Bonche, 

Nucl. Phys. A621, 706 (1997). 
[68] J. Terasaki and J. Engel, Phys. Rev. C 74, 044301 (2006). 
[69] A. Ansari and P. Ring, Phys. Rev. C 74, 054313 (2006). 
[70] A. P. Severyukhin, V. V. Voronov, and N. Van Giai, 

Phys. Rev. C 77, 024322 (2008). 
[71] A. P. Severyukhin, V. V. Voronov, and N. Van Giai, Eur. 

Phys. J. A 22, 397 (2004). 
[72] S. Raman, C. W. Nestor, and P. Tikkanen At. Data Nucl. 

Data Tables 78, 1 (2001). 
[73] D. C. Radford et ai, Nucl. Phys. A752, 264c (2005). 
[74] A. Banu et al, Phys. Rev. C 72, 061305 (2005). 
[75] C. Vaman et ai, Phys. Rev. Lett. 99, 162501 (2007). 
[76] J. Cederkall et al, Phys. Rev. Lett. 98, 172501 (2007). 
[77] A. Ekstrom et ai, Phys. Rev. Lett. 101, 012502 (2008). 
[78] M. A. Kennedy, P. D. Cottle, and K. W. Kemper, Phys. 

Rev. C 46, 1811 (1992). 



